----------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  Muris-FEOL-Illustration.log
  log type:  text
 opened on:   26 Nov 2015, 16:59:08

{ [code omitted] }



. 
. qui eststo: xtologit RTHLTH FAMINCY AGEINC FAMSY

. lincom [cut2]_cons-[cut1]_cons

 ( 1)  - [cut1]_cons + [cut2]_cons = 0

------------------------------------------------------------------------------
      RTHLTH |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         (1) |   2.010345   .0451644    44.51   0.000     1.921825    2.098866
------------------------------------------------------------------------------

. lincom [cut3]_cons-[cut2]_cons

 ( 1)  - [cut2]_cons + [cut3]_cons = 0

------------------------------------------------------------------------------
      RTHLTH |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         (1) |   2.955451   .0860579    34.34   0.000     2.786781    3.124122
------------------------------------------------------------------------------

. lincom [cut4]_cons-[cut3]_cons

 ( 1)  - [cut3]_cons + [cut4]_cons = 0

------------------------------------------------------------------------------
      RTHLTH |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         (1) |   2.735264    .230982    11.84   0.000     2.282547     3.18798
------------------------------------------------------------------------------

. 
. qui eststo: xtologit RTHLTH FAMINCY AGEINC FAMSY AGE SEX MSA i.REGION RACEV1X FAMINCY_AV AGEINC_AV FAMSY_AV

. lincom [cut2]_cons-[cut1]_cons

 ( 1)  - [cut1]_cons + [cut2]_cons = 0

------------------------------------------------------------------------------
      RTHLTH |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         (1) |   2.020033   .0453529    44.54   0.000     1.931143    2.108923
------------------------------------------------------------------------------

. lincom [cut3]_cons-[cut2]_cons

 ( 1)  - [cut2]_cons + [cut3]_cons = 0

------------------------------------------------------------------------------
      RTHLTH |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         (1) |   2.970314   .0864621    34.35   0.000     2.800851    3.139777
------------------------------------------------------------------------------

. lincom [cut4]_cons-[cut3]_cons

 ( 1)  - [cut3]_cons + [cut4]_cons = 0

------------------------------------------------------------------------------
      RTHLTH |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         (1) |   2.740316   .2309816    11.86   0.000       2.2876    3.193031
------------------------------------------------------------------------------


{ [code omitted] }
. 
. // Present
. sum prob* if RTHLTH==2 & AGEOLD==15

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
       prob1 |       133    .4800521       .1328   .1107167   .7555807
       prob2 |       133    .3784125    .0622417   .2032742    .466035
       prob3 |       133    .1325365    .0738455   .0389492   .4640335
       prob4 |       133    .0084099    .0066344   .0020539   .0482983
       prob5 |       133     .000589    .0004767    .000142   .0035149
-------------+--------------------------------------------------------
    probnew1 |       133     .473685    .1323986   .1080434    .750477
    probnew2 |       133    .3815947    .0610447   .2072814   .4660232
    probnew3 |       133    .1354748    .0749066   .0399847    .469519
    probnew4 |       133    .0086402     .006811   .0021109   .0495658
    probnew5 |       133    .0006054      .00049    .000146   .0036123


{ [code omitted] }

. //
. // Column 3
. // Estimator in the paper uses all transformations
. //
. eststo: clogit ybin `RHS', group(newid) cluster(CASEID)
note: multiple positive outcomes within groups encountered.
note: 49046 groups (98044 obs) dropped because of all positive or
      all negative outcomes.

Iteration 0:   log pseudolikelihood = -3554.0933  
Iteration 1:   log pseudolikelihood =  -3534.318  
Iteration 2:   log pseudolikelihood = -3534.2244  
Iteration 3:   log pseudolikelihood = -3534.2244  

Conditional (fixed-effects) logistic regression   Number of obs   =      34196
                                                  Wald chi2(6)    =    2010.99
                                                  Prob > chi2     =     0.0000
Log pseudolikelihood = -3534.2244                 Pseudo R2       =     0.7018

                              (Std. Err. adjusted for 2689 clusters in CASEID)
------------------------------------------------------------------------------
             |               Robust
        ybin |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     FAMINCY |   .0635977   .0758006     0.84   0.401    -.0849688    .2121642
      AGEINC |  -.0351112   .0162837    -2.16   0.031    -.0670267   -.0031958
       FAMSY |   .2285883   .1583783     1.44   0.149    -.0818275     .539004
             |
         cut |
          2  |   1.873126   .0464733    40.31   0.000      1.78204    1.964212
          3  |   4.789859   .1327821    36.07   0.000     4.529611    5.050107
          4  |    7.40066   .3297054    22.45   0.000      6.75445    8.046871
------------------------------------------------------------------------------
(est5 stored)

. lincom 3.cut-2.cut

 ( 1)  - [ybin]2.cut + [ybin]3.cut = 0

------------------------------------------------------------------------------
        ybin |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         (1) |   2.916733   .1168164    24.97   0.000     2.687777    3.145689
------------------------------------------------------------------------------

. lincom 4.cut-3.cut

 ( 1)  - [ybin]3.cut + [ybin]4.cut = 0

------------------------------------------------------------------------------
        ybin |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         (1) |   2.610802    .284268     9.18   0.000     2.053647    3.167957
------------------------------------------------------------------------------

. 
. 
. ///////////////////////////////////
> // 5. Present the results  ////////
. ///////////////////////////////////
> 
. di "Remember: reported results for transformed estimators should be multiplied by (-1)."
Remember: reported results for transformed estimators should be multiplied by (-1).

. 
. esttab, se b(a2) mtitle("pi=(3,3)" "RE" "CRE" "BUC" "CLE")

--------------------------------------------------------------------------------------------
                      (1)             (2)             (3)             (4)             (5)   
                 pi=(3,3)              RE             CRE             BUC             CLE   
--------------------------------------------------------------------------------------------
main                                                                                        
FAMINCY             -0.33           -0.38***        -0.10           0.090           0.064   
                   (0.26)         (0.034)         (0.062)         (0.070)         (0.076)   

AGEINC             -0.083          -0.014           0.017          -0.021          -0.035*  
                  (0.049)        (0.0072)         (0.013)         (0.015)         (0.016)   

FAMSY                0.46           0.094*          -0.19            0.20            0.23   
                   (0.60)         (0.040)          (0.14)          (0.15)          (0.16)   

AGE                                                 0.043***                                
                                                 (0.0082)                                   

SEX                                                 -0.23**                                 
                                                  (0.074)                                   

MSA                                                 -0.11                                   
                                                   (0.11)                                   

1.REGION                                                0                                   
                                                      (.)                                   

2.REGION                                           -0.057                                   
                                                   (0.13)                                   

3.REGION                                           -0.067                                   
                                                   (0.11)                                   

4.REGION                                            -0.12                                   
                                                   (0.12)                                   

RACEV1X                                           -0.0100                                   
                                                  (0.034)                                   

FAMINCY_AV                                          -0.41***                                
                                                  (0.075)                                   

AGEINC_AV                                          -0.041**                                 
                                                  (0.016)                                   

FAMSY_AV                                             0.31*                                  
                                                   (0.15)                                   

1.cut                                                                   0               0   
                                                                      (.)             (.)   

2.cut                                                                   0            1.87***
                                                                      (.)         (0.046)   

3.cut                                                                   0            4.79***
                                                                      (.)          (0.13)   

4.cut                                                                   0            7.40***
                                                                      (.)          (0.33)   
--------------------------------------------------------------------------------------------
cut1                                                                                        
_cons                               0.037           -0.49**                                 
                                  (0.040)          (0.18)                                   
--------------------------------------------------------------------------------------------
cut2                                                                                        
_cons                                2.05***         1.53***                                
                                  (0.053)          (0.19)                                   
--------------------------------------------------------------------------------------------
cut3                                                                                        
_cons                                5.00***         4.50***                                
                                   (0.11)          (0.20)                                   
--------------------------------------------------------------------------------------------
cut4                                                                                        
_cons                                7.74***         7.24***                                
                                   (0.26)          (0.31)                                   
--------------------------------------------------------------------------------------------
sigma2_u                                                                                    
_cons                                3.11***         3.09***                                
                                   (0.20)          (0.20)                                   
--------------------------------------------------------------------------------------------
N                     318            8265            8265            4730           34196   
--------------------------------------------------------------------------------------------
Standard errors in parentheses
* p<0.05, ** p<0.01, *** p<0.001

. eststo clear

. 
. log close
      name:  <unnamed>
       log:  Muris-FEOL-Illustration.log
  log type:  text
 closed on:   26 Nov 2015, 16:59:42
----------------------------------------------------------------------------------------------------------------------------------------------
